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Abstract 

We consider the dynamics of rotational excitations placed on a single molecule in spatially dis¬ 
ordered ID, 2D and 3D ensembles of ultracold molecules trapped in optical lattices. The disorder 
arises from incomplete populations of optical lattices with molecules. This leads to a model cor¬ 
responding to a quantum particle with long-range tunnelling amplitudes moving on a lattice with 
the same on-site energy but with forbidden access to random sites (vacancies). We examine the 
time and length scales of Anderson localization for this type of disorder with realistic experimen¬ 
tal parameters in the Hamiltonian. We show that for an experimentally realized system of KRb 
molecules on an optical lattice this type of disorder leads to disorder-induced localization in ID 
and 2D systems on a time scale t ~ 1 sec. For 3D lattices with 55 sites in each dimension and 
vacancy concentration 90 %, the rotational excitations diffuse to the edges of the lattice and show 
no signature of Anderson localization. We examine the role of the long-range tunnelling ampli¬ 
tudes allowing for transfer of rotational excitations between distant lattice sites. Our results show 
that the long-range tunnelling has little impact on the dynamics in the diffusive regime but affects 
significantly the localization dynamics in lattices with large concentrations of vacancies, enhancing 
the width of the localized distributions in 2D lattices by more than a factor of 2. Our results raise 
a general question whether quantum particles with long-range tunnelling can undergo quantum 
localization in 3D lattices with substitutional disorder. 
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I. INTRODUCTION 


Scattering of electrons by impurities in disordered crystals leads to Anderson localization 
[I; £2], which determines the conductor - insulator transitions in metals [3], quasi-crystals |J] 
and granular metal films [5] , and is associated with many interesting phenomena, such as the 
quantized phases of the integer Hall effect |2UBj- If written in a second-quantized form, the 
Hamiltonian of an electron in a disordered lattice can be mapped onto the Hamiltonians for a 
variety of systems [7]. This mapping can be used to study transport properties of electrons 
in disordered lattices by examining the properties of other quantum particles and quasi¬ 
particles. For example, disordered-induced localization has been studied with microwaves 
in a tubular waveguide m, optical photons in opaque media mm and ultracold atoms 
in optical lattices mm- Despite these studies, there are still many open questions about 
Anderson localization. For example, the scaling hypothesis providing the relation between 
quantum localization in systems of different dimensionality |16] still awaits experimental 
confirmation. Equally important and widely debated are the effects of inter-particle inter¬ 
actions nn and dissipative forces [T8] on localization in far-from-equilibrium systems or 
the role of long-range tunnelling amplitudes P2H2D allowing particles to transition between 
distant lattice sites. 

The development of experimental methods for producing molecules at ultracold temper¬ 
atures has created new possibilities for studying disorder-induced localization. As demon¬ 
strated in recent experiments [22H2S], an ensemble of ultracold KRb diatoms, produced by 
photoassociation of ultracold K and Rb atoms and trapped in an optical lattice, forms a ran¬ 
dom spatial distribution of KRb molecules in the ro-vibrational ground state. If the trapping 
field of the optical lattice has high intensity, the molecular motion is strongly confined and 
the spatial distribution is quasi-stationary, eliminating molecule - molecule collisions. The 
excitation of molecules from the rotational ground state to the lowest rotationally excited 
state leads to the formation of delocalized excitations [23 [2Ej, representing the wave packets 
of rotational Frenkel excitons mm- The purpose of the present work is to explore if these 
wave packets can be used as probe particles to study dynamics of localization induced by 
the disorder potential stemming from the randomness of the molecular positions. 

The rotational excitations in an ensemble of polar molecules randomly distributed on an 
optical lattice possess several unique properties: 
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• The rotationally excited states of molecules have long radiative lifetimes (> 1 — 10 
sec), allowing one to probe coherent dynamics over long time scales. 

• The rotational excitations can be placed on molecules in well defined lattice sites. This 
can be achieved by applying a spatial gradient of an electric field and a microwave 
field pulse resonant with rotational transitions for molecules in the desired part of the 
lattice pl|. The rotational excitations thus produced form spatially localized wave 
packets [28]. Once the gradient of the field is removed, these wave packets can travel 
throughout lattice due to resonant energy transfer between molecules, undergoing 
quantum walk (HUSH and scattering by empty lattice sites. By shaping the gradient 
of the field, it may be possible to prepare the excitation wave packets representing 
superpositions of excitations in spatially separated lattice sites, thus enabling the 
controlled study of the role of initial conditions and entanglement between spatially 
separated particles on the dynamics of quantum walk. 

• The interactions between the rotational excitations can be tuned by applying a DC 
electric field [ 32] , 

• The transfer of rotational excitations between molecules is enabled by long-range dipole 
- dipole interactions, leading to long-range transfer of excitations between distant 
lattice sites. This can be used to study the effects of long-range tunnelling on disorder- 
induced localization. 

• The translational motion of molecules can be controlled by varying the strength of 
the optical lattice potential. As shown in Ref. [33] 133] . the translational motion of 
molecules can be coupled to the rotational excitations. This can be used to study 
disorder-induced localization and diffusion of quantum particles in the presence of 
controlled coupling to a dissipative environment. 

• The optical lattice potential can be designed to produce a three-dimensional (3D) 
finite lattice of molecules or to separate molecules in one or two dimensions, effectively 
confining the rotational excitations to one-dimensional (ID) or two-dimensional (2D) 
disordered lattices [33]. 

These properties make the rotational excitations in spatially disordered ensembles of ultra¬ 
cold molecules a unique new platform for the study of quantum walk and disorder-induced 
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localization. In particular, the ability to change the dimensionality of the optical lattice 
potential opens up a unique possibility to verify the universality of Anderson localization in 
different dimensionalities [T6] and the ability to tune the interactions between the rotational 
excitations ra suggests a new way to study the effects of particle interactions on Anderson 
localization. However, the microscopic Hamiltonian describing rotational excitations in an 
optical lattice leads to a specific model corresponding to a quantum particle moving on a 
lattice with the same on-site energy but with forbidden access to random sites, with the 
tunnelling amplitudes decaying as the inverse cube of the hopping distance. Whether or not 
such a model can be used to observe Anderson localization in finite, or even infinite [SUEZ!, 
2D systems and whether or not such a model allows Anderson localization in 3D systems of 
any size [38] is - to the best of our knowledge - unknown [39] . 

The specific goal of this work is to consider the following questions: 

(i) While quantum particles placed in a 2D disordered potential are generally expected 
to undergo Anderson localization heissiezi, the localization length is sensitive to 
the microscopic details of the particle - disorder interactions and can be very large. 
Given that the experiments with optical lattices employ finite ensembles of molecules, 
typically with less than 100 sites per dimension, is the disorder potential for rota¬ 
tional excitations in a partially populated optical lattice strong enough to allow for 
observations of Anderson localization in 2D systems on the length scale of 100 lattice 
sites? 

(ii) If the disorder potential is strong enough to induce Anderson localization in finite 2D 
arrays, what is the minimum concentration of lattice vacancies leading to localization 
lengths less than 100 lattice sites? 

(iii) Quantum particles placed in a 3D disordered potential may or may not undergo An¬ 
derson localization pT6| , depending on the microscopic details of the particle - disorder 
interactions. Is the disorder potential for rotational excitations in 3D disordered en¬ 
sembles strong enough to induce Anderson localization? 

(iv) What is the effect of the long-range tunnelling amplitudes on the dynamics in the 
disordered system of this specific kind? 
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(v) What is the time scale of the disorder-induced localization, if any, for excitations 
placed in specific lattice sites? 

If the disorder potential is too weak, the localization dynamics may not be observ¬ 
able over experimentally feasible time scales. If the disorder potential is too strong, 
localized rotational wave packets may remain pinned to their original location and not 
exhibit any dynamics of quantum walk over experimentally feasible time scales. 

Questions (i) - (iv) are general for the lattice model with tunnelling amplitudes randomly 
omitted and the results pertaining to question (v) can be renormalized for a particular 
system by re-scaling the time dependence of observables. To answer these questions, we 
consider the dynamics of rotational excitations placed in individual sites of an optical lat¬ 
tice partially populated with KRb molecules. To simulate an experiment with destructive 
measurements, we average the dynamics of rotational energy transfer between molecules in 
the disordered ensembles over multiple disorder orientations and examine the time scales of 
the formation of disorder-induced distributions of rotational excitations in lattices of various 
dimensionality. We examine the effects of the long-range tunnelling amplitudes in lattices of 
various dimensionality and discuss the effect of disorder on diffusion of rotational excitations 
in lattices of various dimensionality. 

II. METHODOLOGY 

The Hamiltonian for an ensemble of molecules on an optical lattice of arbitrary dimen¬ 
sionality can be written as 


n n n' 

where H n is the Hamiltonian of an isolated molecule placed in lattice site n positioned 
at r = na, a is the lattice constant, and V nn ‘ is the dipole - dipole interaction between 
molecules in sites n and n' . The vector index is n = n x for ID lattices, n = (n x ,n y ) for 2D 
lattices and n = (n x ,n y ,n z ) for 3D lattices, with each of n x , n y and n z running from — N 
to N, so that J\f = 2 N + 1 is the number of lattice sites per dimension. We assume that 
the translational motion of molecules can be neglected. This approximation is valid for the 
trapping held potential with the trapping frequency exceeding 100 kHz |33j . 
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Following the experimental work in Ref. [25], we assume that each molecule is initially 
in the lowest eigenstate of H n , corresponding to zero rotational angular momentum j — 0 of 
the molecule. Some of the molecules are subsequently excited to an isolated hyperhne state 
within the triplet of the hyperhne states corresponding to the rotational angular momentum 
j = 1 and the projection of j on the space-fixed quantization axis rrij = —1. The lowest 
energy state and this specific excited state for the molecule in site n are hereafter denoted 
as | g n ) and \e n ). The states | g) and \e) have the opposite parity. 

Since the dipole - dipole interaction operator can be written as a sum over products 
of rank-1 spherical tensors acting on the subspace of the individual molecules, the matrix 
elements {g n \{g n '\Vn,n'\g n )\gn') vanish so the molecules are non-interacting when in the 
ground state. However, the dipole - dipole interaction has non-zero matrix elements between 
the product states \g n )\e n /) and \e n )\g n i). The matrix elements 

tn,n' {gn | (fin' | ^n,n' | ) (2) 


stimulate resonant energy transfer of the | g) —y |e) excitation between molecules in sites n 
and n' . For the specific states | g) = | j = 0) and |e) = | j = 1, rrij = —1), the matrix elements 
(|2| have the following form: 


tn,n' 


(Iri/lr 


47reo 6a 3 1 n — n'| 3 


(3 cos 2 9 nn > - 1) 


( 3 ) 


where d n = d ri > = 0.57 Debye is the permanent dipole moment of KRb molecules ESI j 9 nn i is 
the angle between the intermolecular axis joining the molecules in lattice sites n and n' and 
the ^-axis. We assume that 6 nn / = 7t/2 for ID and 2D lattices. The molecules are trapped 
in a lattice with lattice constant a = 532 nm [25]. For the specific system considered here, 
the value of these matrix elements is 




(3 COS^ 0 rin* 1) 


\n — n 


/ 13 


( 4 ) 


with a = 52.12 Hz. 

The matrix elements (|2]) also lead to the delocalized character of the rotational excitations. 


6 





In general, the | g) —>• |e) excitation generates the following many-body excited state: 

IVw) = \e n ) JJ I gi). ( 5 ) 

n i^n 

In an ideal lattice, the coefficients C n = e} ap ' n j\[M and l^w) => h/W,(p)) represents a 
rotational Frenkel exciton with wave vector p [5DJ , i.e. a rotational excitation completely 
delocalized over the entire lattice. If the lattice is disordered, the excited state is a localized 
coherent superposition of the Frenkel excitons with different p [28]. 

When the lattice is partially filled with molecules, the excitation energy of the molecular 
ensemble depends on the location of molecules in the excited state. The excited states of the 
many-body system of molecules in a disordered lattice are the eigenstates of the following 
Hamiltonian: 


H ^ T ^ ^ ^ ^ t. n n ', (6) 

n n n' 

where the operator c n removes an excitation from site n, w n is the excitation energy of a 
molecule in site n, if the site is populated, and zero otherwise. The transition amplitudes 
tn,n' are given by Eq. (|2]) if both of the sites n and n' are populated by molecules and zero 
otherwise. It is this disorder in the transition amplitudes t n n i that localizes the delocalized 
excitations. 

In order to compute the time-evolution of the rotational excitations, we diagonalize the 
matrix of the Hamiltonian (Jb]) in the site-representation basis |n) = cjj 0), where the vacuum 
state is 


io>=n is-). w 

m 

with both m and n running only over the lattice sites populated with molecules. 

The eigenstates of the full Hamiltonian (J6]) , 


n 


( 8 ) 
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are then used to compute the time evolution of the excitation wave packet as follows 


\m) = E Cxe- iE ^ l \X) = £/»(«) |n), 


( 9 ) 


A 


n 


where E\ are the eigenvalues of the Hamiltonian (J6]), the coefficients 


/„(«) = 


( 10 ) 


A 


and the coefficients C\ are the projections of the excitation wave packet at t — 0 onto the 
eigenstates of the Hamiltonian, 


c x = (\m = o)). 


(ii) 


The quantities |/ n (t)| 2 are averaged over multiple calculations with different realizations of 
disorder (i.e. different random distributions of empty lattice sites). 

As will be demonstrated in the following section, the rotational excitations scattered by 
the empty lattice sites form peaked distributions. In this article, we will use two measures 
of distributions: the distribution width L and the standard deviation cr r . The distribution 
width L is defined as the length of the lattice (for a ID lattice), the radius of a circle (for 
a 2D lattice) or the radius of a sphere (for a 3D lattice) containing 90 % of the excitation 
probability. The standard deviation is defined as ay = yj (r 2 ) — (r) 2 , where r 2 = x 2 for 
a ID lattice, r 2 = x 2 + y 2 for a 2D lattice and r 2 = x 2 + y 2 + z 2 for a 3D lattice. The 
distribution width characterizes the physical spread of the rotational excitation over the 
molecular ensemble. However, this quantity may be significantly affected by fluctuations 
due to specific disorder realizations. The standard deviation depends on the shape of the 
distribution but is less affected by the fluctuations. We will use the standard deviation in 
order to characterize non-ambiguously the dynamical properties such as the time dependence 
of the distributions. 

III. QUANTUM WALK OF ROTATIONAL EXCITATIONS 

As originally proposed by DeMillc [59], ultracold molecules trapped in an optical lattice 
can be selectively excited by applying a gradient of a DC electric held that separates the 
rotational energy levels of molecules in different lattice sites by a different magnitude. This 



method can, in principle, be used to excite rotationally a single molecule in a specific lattice 
site. If the gradient of the field is removed, the rotational excitation thus generated may 
be transferred to molecules in other lattice sites by resonant energy transfer. The dynamics 
of this energy transfer can be probed by applying another gradient of the DC electric field 
at time t and detecting either the excited molecules or molecules remaining in the ground 
state using resonantly-enhanced multiphoton ionization (REMPI). The REMPI detection 
is a destructive measurement. If the optical lattice is partially populated with randomly 
distributed molecules, repeating the experiment multiple times to determine the dependence 
of rotational energy transfer on t is equivalent to averaging over different realizations of 
disorder (i.e. different distributions of empty lattice sites). 

In order to simulate the outcome of such measurements, we consider an optical lattice 
partially populated with randomly distributed KRb molecules. We compute the time evolu¬ 
tion of a rotational excitation placed in a single lattice site at t = 0 and average the results 
at each t over multiple random distributions of empty lattice sites with a fixed concentra¬ 
tion of vacancies. Figures 1, 2 and 3 illustrate the characteristic distributions of rotational 
excitations as a result of quantum walk of a single rotational excitation initially placed in 
the centre of a ID, 2D (square) or 3D (cubic) lattice, respectively. In the following sec¬ 
tions, we explore the detailed dynamics of how these distributions are formed and the role 
of long-range tunnelling effects in determining the quantum walk dynamics. 

A. Localization dynamics in ID lattice 

A quantum particle placed in a ID disordered lattice must undergo Anderson localization 
[T]. Mapped onto a rotational excitation travelling in a disordered ensemble of molecules, 
this well-known result is reflected in the formation of the exponentially localized distribu¬ 
tions shown in Figure 1. These distributions characterize the probability of molecules in 
the corresponding lattice sites to be in the rotationally excited state. Figure 1 illustrates 
that rotational excitations in an optical lattice partially populated with molecules behave as 
quantum particles in a disordered lattice. It is important to note that the distributions pre¬ 
sented in Figure 1 are time-independent. Although the rotational excitation initially placed 
in a specific lattice site represents a wave packet that, given sufficient time, can potentially 
explore the entire lattice and exhibit revivals, the averaging over disorder realizations leads 
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to stationary distributions in the limit of long time. In order to illustrate this, we present 
in Figure 4 the distributions of the rotational excitation computed at different times. 

The lower panel of Figure 4 illustrates the time dependence of the distribution width L 
in a ID lattice with the vacancy concentrations 10 % and 20 %. The results show that for 
the parameters representing rotational excitations in an ensemble of KRb molecules, the 
averaged excitation probability distributions change most significantly at t between zero 
and 100 - 300 ms. The distribution widths exhibit a temporary peak at short times. This 
oscillation survives the averaging over disorder realizations and is present in the results for 
different concentrations of vacancies. Since the distribution width presented in Figure 4 is 
much smaller than the lattice size (1001 x a), this short-time oscillation cannot be due to 
reflection from the lattice boundaries and is likely the result of constructive interference due 
to back-scattering of the rotational excitation from the empty lattice sites. 

As indicated by the results in the lower panel of Figure 4, the approach to the stationary 
distributions is determined by the strength of the disorder potential, i.e. the concentration 
of lattice vacancies. This is illustrated more clearly in the upper panel of Figure 5 showing 
the time dependence of the standard deviations oy for different concentrations of vacancies. 
Figure 5 demonstrates that the rotational excitations in a ID lattice of KRb molecules 
reach the time-independent distributions in less than one second, provided the vacancy 
concentration is > 3%. 

It is instructive to examine the dependence of the minimum time required to form the sta¬ 
tionary distributions on the concentration of vacancies. Figure 6 shows the time r required 
for a rotational excitation initially placed in the center of a lattice to form a time-independent 
distribution as a function of vacancy concentration in lattices of different dimensionality. 
Interestingly, Figure 6 reveals that the dependence of r on the vacancy concentration is 
qualitatively different for lattices of different dimensionality. In ID lattices, r is a monoton- 
ically decaying function of the vacancy concentration. This illustrates that increasing the 
disorder of ID lattices impedes the dynamics of quantum walk leading to narrower distribu¬ 
tions of the rotational excitation (cf., Figure 1) and decreases the time required to reach the 
stationary distributions (less time since the quantum walker has less space to explore). The 
calculations shown in Figure 6 have two implications: (i) the time evolution of the rotational 
excitations in ID disordered lattices can be used as a probe of the disorder strength; (ii) 
the mechanisms for the formation of the distributions shown in Figures 1-3 are different 
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for lattices of different dimensionality. In the following section, we correlate the rising or 
falling dependence of r on the concentration of vacancies with the regime of classical diffu¬ 
sion or sub-diffusion, arguing that the sign of the gradient of r can be used as a signature 
of sub-diffusion stimulated by Anderson localization. 

B. Localization and diffusion in 2D and 3D lattices 

As can be seen from Figures 2 and 3 and the comparison of the different curves in Figures 
5 and 6, the dynamics of rotational energy transfer between molecules trapped in 2D and 
3D lattices is quantitatively and qualitatively different from the localization dynamics in ID 
lattices. It is apparent that the formation of localized distributions requires higher concen¬ 
trations of vacancies and the approach to time-independent distributions takes longer time 
in lattices of higher dimensionality. In particular, Figure 6 reveals that the dependence of 
the time r required to reach a stationary distribution on the vacancy concentration is quali¬ 
tatively different for lattices of different dimensionality. In ID lattices, r is a monotonically 
decaying function of the vacancy concentration. In 2D lattices, r exhibits a maximum as a 
function of the vacancy concentration. In 3D lattices, r increases monotonically with the 
disorder strength. 

In order to elucidate the dynamics of quantum walk in higher dimensions and the impli¬ 
cations of Figure 6, we present in Figure 7 the time dependence of the distribution widths 
L computed for a rotational excitation in 2D and 3D lattices of varying size with varying 
concentrations of vacant sites. The results for the 2D lattices illustrate that at low con¬ 
centration of vacancies (20 %) corresponding to the range, where the dependence of r on 
the disorder strength rises, the rotational excitation width approaches the size of the lattice 
in the limit of long time. This shows that the localization length of the excitation for this 
concentration of vacancies is greater than the size of the lattice. When the concentration 
of empty sites is 70 %, the rotational excitation in a 2D lattice is well localized and the 
different curves shown in Figure 7 converge to the same value in the limit of long time. This 
localization regime corresponds to the range in Figure 6, where r decreases with the disorder 
strength. 

Since the experiments with ultracold molecules on optical lattices employ finite samples 
of molecules, typically ~ 100 per dimension ra. it is important to find the threshold 
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concentration of vacancies that leads to localized states in 2D lattices with dimension < 
100 x 100. By repeating the calculation of Figure 7 for different concentrations of vacancies, 
we find that L(t —> oo) < 100, when the concentration of vacancies exceeds 50%. This is 
consistent with the results in Figure 6, which shows that the dependence of r on the vacancy 
concentration exhibits a maximum at the concentration of vacancies between 40% and 50%. 
We conclude that the rotational excitations can be localized in 2D disordered lattices of 
molecules with 100 x 100 lattice sites, provided the concentration of vacancies is > 50%. 

The calculations for the rotational excitation dynamics in 3D lattices (lower panel of 
Figure 7) show that, for each of the lattices considered, the distribution width approaches 
the size of the lattice in the limit of long time, even for the concentration of empty sites 
as large as 90 %. This indicates that the rotational excitation in 3D lattices with 90 % of 
empty sites is allowed to diffuse to the edges of the lattice. This is also consistent with the 
results in Figure 6 showing that r rises monotonically for 3D lattices. This indicates that the 
rotational excitations do not exhibit Anderson localization on the length scale considered. 
To prove this, we compare the dynamics of the spatial expansion of the rotational excitations 
in 3D disordered lattices with classical diffusion. 

The diffusion of quantum particles in disordered lattices can be examined for an alterna¬ 
tive signature of Anderson localization. A quantum particle with cosine dispersion placed 
in an individual site of an ideal lattice expands ballistically to form a distribution with the 
standard deviation of = (r 2 ) — (r) 2 oc t 2 [mug. In contrast, the diffusion of a classi¬ 
cal particle in Brownian motion is characterized by of oc t. In the presence of Anderson 
localization, the diffusive dynamics must be suppressed. 

The time-dependence of oy can thus be used as a quantitative measure of quantum 
interference, constructive or destructive, on particle propagation in disordered lattices. 

The upper panels of Figure 8 (2D lattices) and Figure 9 (3D lattices) show the time 
dependence of of computed for 2D and 3D lattices with various concentrations of empty 
lattice sites. As expected, we find that of oc t 2 for ideal lattices, both in 2D and 3D. 
However, the time-dependence of of is dramatically modified by the presence of vacancies. 
In particular, Figure 8 illustrates that the rotational excitation in disordered lattices expands 
ballistically at short times with of oc t 2 and exhibits the dependence of oc t, characteristic of 
the diffusive regime, at later times. That of becomes a linear function of time illustrates that 
the presence of vacancies destroys interferences accelerating diffusion of quantum particles 
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[I3til6] . At even later times (t > 300 ms), the dependence of of on time departs from 
linearity and the propagation of the rotational excitations enters a sub-diffusive regime 
(of oc t 7 with 7 < 1 ). This can happen either because of Anderson localization or because 
of the reflection from the boundaries. 

Since the calculations of Figures 8 and 9 are performed for finite lattices, it is important 
to examine the effects of the boundaries. In particular, it is important to determine if the 
departure of of from linearity is due to disorder or due to the finite size effects. In order to 
do this, we first consider the ballistic expansion of the rotational excitations in ideal lattices 
of varying size (upper right panels of Figures 8 and 9). We observe that, for each lattice size, 
of increases as a quadratic function of time at t < t\, until the reflection from the boundaries 
slows down the ballistic expansion. We denote by t\ the time, at which the dependence of 
of on time in an ideal lattice departs from the quadratic function due to the reflection from 
the boundaries. The insets of Figures 8 and 9 show the dependence of t\ on the lattice 
size. In the presence of disorder, the expansion is significantly slower and the reflection from 
the boundaries must occur at times later than t\. Therefore, we can use the time fy as the 
lower limit for the onset of the boundary effects and conclude that the dynamics, whether 
in disordered or ideal lattices, is not affected by the boundary effects at t < t\. 

The lower panel of Figure 8 shows that the dynamics of rotational excitations in 2D 
lattices with the disorder concentration > 50% enters the sub-diffusive regime at times 
t < ti. This proves that the transition to the sub-diffusive regimes at these concentration 
of vacancies is due to disorder-induced localization, and not due to the boundary effects. 
For the disorder concentration 20%, we observe that the departure of of from linearity 
occurs at t > fy. These observations are consistent with the results of Figures 6 and 7. 
We thus conclude that the increase of the time required for the formation of the stationary 
distributions shown in Figure 6 reflects that increasing the disorder strength slows down the 
diffusion of rotational excitations to the lattice boundaries, hence the dramatically different 
dependence of r on the disorder strength for 2D (at low vacancy concentrations) and 3D 
lattices from that in ID lattices. 

For 3D lattices, the transition to the sub-diffusive regimes always occurs at £ > A, 
indicating (but not proving) that the departure of of from linearity is due to the finite size 
of the system, and not the scattering by disorder. The lower panel of Figure 9 shows the 
results for an extreme case of a 3D lattice with 90 % of the lattice sites vacant. To prove 
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that the disorder does not lead to the sub-diffusive regime in 3D lattices, we perform the 
calculations illustrated in the lower panel of Figure 9 for different lattice sizes. For each 
lattice size, we determine the time i 2) at which the expansion of a 2 undergoes the transition 
to the sub-diffusive regime. We find that the dependence of U on the lattice size is a linear 
function of the lattice size for extended 3D lattices (31 3 - 55 3 sites), as expected for a particle 
in the diffusive regime. This proves that there is no Anderson localization in 3D lattices on 
the length scale of 55 3 lattice sites for the disordered system considered here. 

C. Effect of long-range tunnelling 

Most models of quantum particles in disordered lattices are based on the nearest neigh¬ 
bour (tight-binding) approximation assuming that the tunnelling amplitudes t n n i in Eq. ([6]) 
are non-zero only when n and n' are adjacent lattice sites. However, the tunnelling ampli¬ 
tudes (J 2 | responsible for rotational energy transfer between molecules decay as oc l/\n — n'[’ 
functions of the lattice site separation. This makes the rotational excitations in an ensemble 
of polar molecules a new platform for the experimental investigation of long-range tunnelling 
effects on the formation of localized states. The 1/| n — n' | 3 dependence of the tunnelling 
amplitudes is particularly interesting [19]. Beyond allowing for direct transitions between 
distant lattice sites, it leads to non-analyticity of particle dispersion at low particle velocities 
ESj. It may also lead to interesting long-range correlations in the localization dynamics of 
multiple interacting particles [47j 4 B] • 

The effect of the long-range character of the tunnelling amplitudes t n ^ n i can be - in 
principle - investigated by comparing the dynamics of rotational excitations in ensembles 
of polar molecules with those in ensembles of homonuclear diatomic molecules. Ultracold 
homonuclear molecules have been produced in multiple experiments worldwide [49J. The 
rotational excitations in an ensemble of non-polar molecules can be resonantly transferred 
between molecules due to quadrupole - quadrupole interactions, the leading term in the 
multipole expansion of the interaction between non-polar molecules separated by a large 
distance. The quadrupole - quadrupole interaction decays as a oc 1/R 5 function of the 
molecule - molecule separation R. Therefore, the rotational excitations in an ensemble of 
homonuclear molecules trapped in an optical lattice should be expected to undergo quantum 
walks governed by the Hamiltonian (| 6 ]) with f nn / oc l/|n — n'\ 5 . 


14 


In this section, we compare the localization dynamics of particles with the tunnelling 
amplitudes t n n > oc l/\n — n'\ 3 and t n n > oc l/\n — n'\ 5 . To remove ambiguity associated with 
the different magnitudes of the dipole and quadrupole moments, we compare the dynamics 
of rotational excitations stimulated by the tunnelling amplitudes as given in Eq. (J2| and 
the tunnelling amplitudes t n , n > = b/\n — n'\ 5 , where the constant b is chosen such that 
tn,n' = t n ,n' for \n — n'\ = 1. In other words, we enhance the quadrupole - quadrupole 
interaction to be the same as the dipole - dipole interaction between molecules in adjacent 
lattice sites so that the difference between the calculations with t n n i and t n n > arises solely 
from the different radial dependence of the tunnelling amplitudes. 

The results presented in Figure 10 illustrate that the long-range matrix elements t n ri j lead 
to less localized distributions in all dimensions. The effect of the long-range matrix elements 
appears to depend both on the strength of the disorder potential and the dimensionality. 
To elucidate these dependencies, we present in Figure 11 the ratio of the distribution widths 
computed in the limit of long time with the long-range amplitudes t n n i oc l/|n — n'| 3 and 
the short-range amplitudes t nn i oc l/\n — n'\ 5 . The long-range matrix elements appear to 
have little effect at low concentrations of empty sites, but increase the distribution widths 
by a factor of 2 or more at large concentrations of vacancies. 

The results presented in Figure 11 indicate that the long-range tunnelling amplitudes 
t n ,ri play little role at low disorder strengths in lattices of higher dimensionality, i.e. when 
the dynamics is expected to be more dominated by diffusion. At the same time, Figure 
11 shows that the long-range matrix elements become exceedingly important in the regime 
where disorder-induced localization is expected to dominate. In order to illustrate this more 
clearly, we show in Figure 12 the distribution widths L computed as functions of time for 2D 
lattices with two concentrations of vacancies: 20 % corresponding to the regime dominated 
by diffusion (cf. Figures 6 and 8) and 70 % corresponding to the regime, where the rotational 
excitations in 2D lattices are well localized. The results show that the long-range tunnelling 
amplitudes insignificantly accelerate the expansion dynamics of the rotational excitation in 
the diffusive regime, enhancing the distribution widths at short times by about 20 % or less 
(upper panel of Figure 12). In stark contrast, the distribution widths obtained with the two 
types of tunnelling amplitudes are dramatically different in the regime of strong localization 
(lower panel of Figure 12). 
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IV. CONCLUSION 


In this work, we study the dynamics of quantum walk of rotational excitations in finite 
disordered ID, 2D and 3D ensembles of ultracold molecules trapped in optical lattices. The 
disorder arises from incomplete populations of optical lattices with molecules. Thus, the 
rotational excitations travelling between molecules in a partially populated optical lattice 
are described by the Hamiltonian (|6]) , which can be obtained from the tight-binding lattice 
model by randomly omitting some of the tunnelling amplitudes and extending the remaining 
tunnelling amplitudes to decay as oc a/\n—n'\ 3 . While the dynamics of particles described by 
such a Hamiltonian should be expected to exhibit the general features of quantum particles 
in a disordered potential, whether or not such a model can be used to observe Anderson 
localization in finite 2D systems and whether or not such a model can be used to observe 
Anderson localization in 3D systems (finite or infinite), depends on the microscopic details 
of particle - disorder interactions. We note that there are only two free parameters in this 
Hamiltonian: the magnitude of the constant a in Eq. Q and the concentration of empty 
lattice sites. While the magnitude of a determines the time scales of quantum walk and 
localization dynamics, the localization length in ID and 2D systems and the presence or 
absence of Anderson localization in 3D systems is entirely determined by the concentration 
of vacancies. 

In order to guide the ongoing experiments with ultracold molecules, we chose the value of 
a to correspond to a specific experimentally realized system of KRb molecules on an optical 
lattice with lattice constant a = 532 nm. We have confirmed that increasing the value of a 
accelerates the time evolution of the dynamics but leads to the same results in the limit of 
long time. The qualitative features of the time-dependent results and the results in the limit 
of long time are, therefore, general for any system described by the Hamiltonian considered 
here. The time-dependent results can be renormalized for other systems by changing the 
value of the coupling constant in Eq. (J2]), which controls the time-evolution of the system. 

The main results of this work can be summarized as following: 

• For the specific system of ID arrays of KRb molecules on an optical lattice with lattice 
constant a = 532 nm, the rotational excitations placed in individual lattice sites form 
localized distributions within t ~ 1 sec. With only 10 % of vacancies, the width of the 
exponentially localized distributions is < 40 a. This indicates that Anderson localiza- 
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tion can be studied with rotational excitations in ID arrays of molecules, providing 
coherence can be preserved for longer than 1 second. 

• In 2D, this particular type of disorder requires concentrations of empty lattice sites 
> 50% in order to localize the rotational excitations within 100 x 100 sites. 

• For 3D disordered arrays with 55 sites in each dimension, we observe no localization 
at vacancy concentrations < 90 %, i.e. the rotational excitations generated in the 
middle of the lattice diffuse to the edges of the lattice as classical particles. The 
vacancy concentration 90 % is near the percolation threshold for a 3D network of sites 
ra- We note that the results of Ref. [[38] indicate that the model for the particle with 
nearest-neighbour hopping on a lattice with random on-site energies and substitutional 
disorder does allow quantum localization in 3D. Figure 5 of Ref. [38] suggests that 
quantum localization is present in a system with substitutional disorder even in the 
absence of disorder in the on-site energies. Our results thus indicate (but not yet prove) 
that the interval of quantum localization in Figure 5 of Ref. [38] between the regimes of 
diffusion and classical localization may be reduced to zero due to long-range tunnelling. 
This raises a general question whether quantum particles with long-range tunnelling 
can undergo quantum localization in 3D lattices with substitutional disorder. 

• Our results show that the long-range character of the tunnelling amplitudes in the 
Hamiltonian ([6]) has little impact on the dynamics of particles in the diffusive regime 
(i.e. in 2D and 3D lattices with low concentrations of vacancies) but affects significantly 
the localization dynamics in lattices with large concentrations of vacancies, enhancing 
the width of the localized particle distributions in 2D lattices by a factor of 2. 

• Our results show that the diffusive vs localization regime can be identified by measuring 
the dependence of the minimum time r (cf. Figure 6) required for a particle placed in 
a specific lattice site to form a time-independent averaged distribution as a function 
of vacancy concentration. 

• Our calculations illustrate that the diffusion of rotational excitations placed in 2D 
lattices with disordered tunnelling amplitudes has three distinct time scales, which 
can be varied by changing the concentration of empty lattice sites. At short times, 
the dynamics is characterized by the quadratic time dependence of of = (r 2 ) — (r) 2 , 
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characteristic of quantum particles with cosine dispersion. The rise of of is suppressed 
by the disorder at later times. This suppression leads to an extended time interval, 
where of exhibits a linear dependence on time, characteristic of classical diffusion. This 
interval of time increases with decreasing disorder strength and can be as long as ~ 100 
ms for the dynamics of in a 2D lattice with 50 % of empty sites. At still longer times, 
the dependence of of on time is sub-diffusive, i.e. slower than linear. This suggests 
that rotational excitations in 2D disordered ensembles of polar molecules can be used 
to study the crossover from ballistic expansion to classical diffusion, where quantum 
interferences enhancing diffusion are suppressed. At longer times, the dynamics of 
rotational excitations in 2D lattices can be used to study the crossover from classical 
diffusion to a sub-diffusive regime, where quantum interferences leading to localization 
become important. 

In our previous work [33], we showed that the rotational excitations of molecules on 
an optical lattice can be coupled to the translational motion of molecules, leading to the 
formation of polarons with a wide range of tunable properties [33]. We also showed that, 
if the molecules on an optical lattice are subjected to DC electric fields, the rotational 
excitations exhibit strong interactions that may - depending on the electric held strength - 
lead to the formation of quasi-particle bound pairs |32j . In combination with this previous 
work, the present results suggest an exciting research avenue for the study of quantum 
diffusion through disordered lattices in the presence of dissipation and the role of inter¬ 
particle interactions on Anderson localization. 
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FIG. 1: Averaged rotational excitation probability distributions formed at t = 2 s by a single 
rotational excitation placed at t = 0 in site n = 0 of a ID lattice of KRb molecules with randomly 
distributed empty sites. The concentration of vacancies is zero (a), 1 % (b), 10 % (c), 20% (d), 
50% (e) and 70% (f). The results are averaged over 100 realizations of disorder. 
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FIG. 2: Averaged rotational excitation probability distributions formed at t = 5 s by a single 
rotational excitation placed at t = 0 in site n x = 0, n y = 0 of a 2D lattice of KRb with a total 
of 51 x 51 sites containing a random distribution of empty sites. The concentration of vacancies 
is zero (a), 10 % (b), 20 % (c), 50% (d), 70% (e) and 85% (f). The results are averaged over 
100 realizations of disorder. For better visualization, the probability values in the two-dimensional 
plots are magnified by the indicated factor. 
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FIG. 3: Averaged rotational excitation probability distributions formed at t = 5 s by a single 
rotational excitation placed at t = 0 in site n x = 0, n y = 0, n z = 0 of a 3D lattice of KRb molecules 
with a total of 31 x 51 x 31 sites containing a random distribution of empty sites. The concentration 
of vacancies is 50% (left panels), 70% (middle panels) and 85% (right panels). The results are 
averaged over 300 realizations of disorder. The upper panels show the part of the distributions for 
x < 0 and y > 0. For better visualization, the probability values in the three-dimensional plots are 
magnified by a factor of 27 (left), 10 (middle), and 2.5 (right). 
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FIG. 4: Upper panel: Time dependence of the rotational excitation probability distributions 
formed from one rotational excitation placed at t = 0 in the lattice site n = 0 of a ID lattice 
of KRb molecules with 20 % of vacancies. Lower panel: Time dependence of the distribution 
width L (in units of the lattice constant a) defined as the range of the lattice containing 90 % 
of the rotational excitation probability: squares - vacancy concentration 10 %; circles - vacancy 
concentration 20 %. The results at each time are averaged over 1000 random realizations of the 
vacancy distributions. 
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FIG. 5: Time dependence of the standard deviations a r of the rotational excitation probability 
distributions formed from one rotational excitation placed at t = 0 in the centre of a lattice partially 
populated with molecules. Upper panel: ID lattice with 1001 sites; middle panel: 2D square lattice 
with 51 x 51 lattice sites; lower panel: 3D cubic lattice with 31 x 31 x 31 lattice sites. The results 
at each time are averaged over > 100 random realizations of disorder. The different sets of data 
correspond to different concentrations of empty lattice sites. 
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FIG. 6: Minimum time r required for a single rotational excitation to form a time-independent 
averaged distribution as a function of vacancy concentration in ID (circles), 2D (squares) and 3D 
(triangles) disordered lattices as a function of vacancy concentration. The results at each vacancy 
concentration are averaged over > 100 realizations of disorder. 
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FIG. 7: Time dependence of the width L (in units of the lattice constant a) of the averaged 
rotational excitation distributions formed by one rotational excitation placed at t = 0 in the centre 
of 2D and 3D lattices with different size. The upper and the lower panels show that the excitations 
diffuse to the edges of the lattice, while the middle plane illustrates that the excitation is localized 
within the lattice and is not affected by the lattice boundaries. 
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FIG. 8: Time dependence of cr 2 = (r 2 ) — (r) 2 for a rotational excitation initially placed in the 
middle of a 2D lattice. Upper left panel: circles - ideal lattice; squares - lattice with 20 % of 
vacancies; triangles - lattice with 50 % of vacancies. Upper right panel: ideal lattice with size 51 
x 51 (diamonds); 101 x 101 (circles); 151 x 151 (triangles); 201 x 201 (squares). Inset: Lattice 
size dependence of time t\ marking the deviation of the wave packet dynamics from that in an 
ideal infinite lattice due to the boundary effects. Lower panel: 2D lattice with 20 % of vacancies 
(squares); 50 % of vacancies (circles); 70 % of vacancies (triangles). The symbols represent the 
numerical calculations; the full curves are the analytical fits cr 2 = Dt 2 ; the dashed curves are the 
linear fits a 2 = bt + c. The size of the lattice for the results in the upper left panel and the lower 
panel is 101 x 101. 
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FIG. 9: Time dependence of cr 2 = (r 2 ) — (r) 2 for a rotational excitation initially placed in the 
middle of a 3D lattice. Upper left panel: 3D ideal lattice (circles), disordered lattice with 50 % 
of vacancies (squares) and disordered lattice with 85 % of vacancies with size 31x31x31 sites. 
Upper right panel: 3D ideal lattice with size 21 x 21 x 21 (squares); 25 x 25 x 25 (triangles); 31 x 
31 x 31 (circles); 35 x 35 x 35 (diamonds). Inset: Lattice size dependence of time t\ marking the 
deviation of the wave packet dynamics from that in an ideal infinitely lattice due to the boundary 
effects. Lower panel: 3D disordered lattice of size 55 x 55 x 55 with 90 % of vacancies. The symbols 
represent the numerical calculations; the full curves are the analytical fits a 2 = Dt 2 -, the dashed 
curves are the linear fits a 2 = bt + c. 
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FIG. 10: Averaged rotational excitation probability distributions formed at t = 2 s (for ID 
lattices) and t = 5 s (for 2D and 3D lattices) by a single rotational excitation placed at t = 0 in 
the center of the lattices. The solid curves are the results of the calculation with the tunnelling 
amplitudes t n n i as defined in Eq. ([2]). The dashed curves are the result of the calculation with the 
tunnelling amplitudes t n ^ oc l/|n — n'| 5 . The middle and lower panels show the cross sections 
of the 2D and 3D distributions along the x-axis. The vertical dotted lines show the distribution 
widths containing 99 % of the rotational excitation. The dimensions of the lattices are 1001 sites 
for ID; 51 X 51 sites for 2D: and 31 X 31 x 31 sites for 3D. 































FIG. 11: The ratio of the widths of the rotational excitation distributions Lg/L l. The distri¬ 
bution widths Lg are computed with the long-range amplitudes t n ^ n i as defined in Eq. ([2]). The 
distribution widths Lg are computed with the short-range amplitudes t nn i oc l/|n — n'\ 5 . The 
results at each vacancy concentration are averaged over > 100 realizations of disorder. 
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FIG. 12: Time dependence of the width L (in units of the lattice constant a) of the averaged 
rotational excitation distributions formed by one rotational excitation placed at t = 0 in the centre 
of 2D lattices with different concentrations of empty lattice sites: circles - the results computed 
with the long-range amplitudes t n n ’ as defined in Eqs. 0 and Q; squares - the results computed 
with the short-range amplitudes t nn > oc 1 /|n — n '\ 5 . 
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